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>~>. We present the clustering measurement of hard X-ray selected AGN in the 

local Universe. We used a sample of 199 sources spectroscopically confirmed de- 
tected by Swift-BAT in its 15-55 keV all-sky survey. We measured the real space 
projected auto-correlation function and detected a signal significant on projected 
scales lower than 200 Mpc/h. We measured a correlation length of ro=5. 56^0^43 
Q ! Mpc/h and a slope 7=1.64lQ;Qy. We also measured the auto-correlation function 

^ ! of Type I and Type II AGN and found higher correlation length for Type I AGN. 

^ ! We have a marginal evidence of luminosity dependent clustering of AGN, as we 

detected a larger correlation length of luminous AGN than that of low luminos- 
ity sources. The corresponding typical host DM halo masses of Swift-BAT are 
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c/5 ! ~ log(M£)A,/H)~ 12-14 h~^M/M0, depending on the subsample. For the whole 

sample we measured log(M£)M//)~ 13.15 h~^M/M0 which is the typical mass 
of a galaxy group. We estimated that the local AGN population has a typical 

^ ! lifetime tagn ~0.7 Gyr, it is powered by SMBH with mass Mbh ~1-10x 10^ M© 

v^ ! and accreting with very low efficiency, log(e)~-2.0. We also conclude that local 
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AGN host galaxies are typically red-massive galaxies with stellar mass of the 
order 2-80x10^° h~^M0. We compared our results with clustering predictions of 



O : merger-dr,ve„ AGN triggermg models aud found a good agreement. 
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Introduction 



It is now commonly believed that almost all galaxies host a central supermassive Black 
Hole (SMBH). Dynamical evidence show that the mass of the central BHs are closely linked 
to the mass as well as the stellar velocity dispersion of the bulge component of the host galaxy 
(Kormendy & Richstone 1995, Magorrian et al. 1998). This suggests that the formation and 
evolution of the spheroidal component of galaxies and their SMBH are closely connected. 
It is of upmost importance to understand the mechanism of funneling interstellar gas into 
the vicinities of the SMBH, triggering the accretion. Galaxy mergers or tidal interaction 
between close pairs may have played a major role (Hopkins et al. 2007), furthermore some 
mechanism internal to the galaxy, like galaxy disk instability may be important. Clustering 
properties of AGN in various redshifts give an important clue to understand which of these 
mechanisms trigger AGN activities in what stage of the evolution of the universe, through, 
e.g., the mass of the Dark Matter Halos (DMH) in which they reside which is linked to BH 
mass life time and Eddington rate. Measurements of the AGN clustering show that AGN 
are typically hosted by DMH with a mass of the order of log(M)~12.0-13.5 M/Mq (Yang 
et al. 2006; Miyaji et al. 2007; Gilli et al. 2009; Coil et al. 2009; Hickox et al. 2009; 
Krumpe et al. 2010). However, these measurements have been produced by using AGN 
samples obtained by optical and soft X-ray (i.e. 0.5-10 keV) surveys. Optical and soft X-ray 
selections miss a major part of the SMBH accretion. In the optical band the selection of 
AGN is biased by galaxy starlight dilution and by dust absorption. Although luminous soft 
X-ray emission is a signature of the presence of an AGN, absorbed sources can be missed 
with a soft X-ray selection, either because they are intrinsically less luminous (Hasinger et 
al. 2008) or because of the high column density. However, X-ray emission from these sources 
leaks out at higher energies (i.e > 5-10 keV) where the efficiency of instruments mounting 
X-ray focusing optics is low. For this reason hard X-ray selected samples could provide clean 
and unbiased samples of AGN. The Swift-BAT all-sky survey provides a spectroscopically 
complete (100 %) sample of local AGN detected in the 15-55 keV energy band, with an 
unprecedented depth and characterization of the source properties, from redshifts to column 
densities. In this letter we present the first study of clustering of hard X-ray selected AGN 
in the local Universe. Throughout this paper we will assume a A-CDM cosmology with 
r2m=0.3, r2A=0.7, Ho=100/i"^ km s~^ Mpc and (78=0.8. Unless otherwise stated errors are 
quoted at the la level. 
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2. The Sample of Swift BAT hard X-ray selected AGN 

The Burst Alert Telescope (BAT; Barthelmy et al. 2005) on board the Swift satellite 
(Gehrels et al. 2004), represents a major improvement in sensitivity for imaging the hard 
X-ray sky. BAT is a coded mask telescope with a wide field of view (FOV, 120° x 90° partially 
coded) aperture sensitive in the 15-200 keV domain. Thanks to its wide FOV and its pointing 
strategy, BAT monitors continuously up to 80% of the sky every day achieving, after several 
years of survey, deep exposure in the entire sky. Results of the BAT survey (Markwardt et 
al. 2005, Ajello et al. 2008, Tueller et al. 2009) show that BAT reaches a sensitivity of 
~1 mCraqj in 1 Ms of exposure. 

The sample used in this work consists of 199 non-blazar AGN detected by BAT during 
the first three years and precisely between March 2005 and March 2008. This sample is 
part of the one used in Ajello et al. (2009) which comprises all sources detected by BAT at 
high (|b|>15°) Galactic latitude and with a signal-to-noise ratio (S/N) exceeding 5 a. The 
reader is referred to Ajello et al. (2009) for more details about the sample and the detection 
procedure. The fiux limit at each direction in the sky has been computed, following Ajello 
et al. (2008), analyzing the local background around that position. For each source we use 
the redshift already provided in Ajello et al. (2009) and the measurement of the absorbing 
column density as determined from joint XMM-Newton/XRT-BAT fits (Burlon et al, in 
preparation). The redshift distribution of the sample is shown in Fig. [T] together with the 
redshift cone diagram of the survey up to 20000 km/s (z~0.07). 



3. The two-point spatial auto-correlation function 

The two-point auto-correlation function (^(r), ACF) describes the excess probability 
over random of finding a pair with an object in the volume dVi and another in the volume dV2, 
separated by a distance r so that dP = n'^[l + C,{r)]dVidV2, where n is the mean space density. 
A known effect when measuring pairs separations is that the peculiar velocities combined with 
the Hubble flow may cause a biased estimate of the distance when using the spectroscopic 
redshift. To avoid this effect we computed the projected ACF (Davies & Peebles 1983): 
w{rp) = 2 Jq""""^ ^('^p; 7i)d7i. Where Vp is the distance component perpendicular to the line of 
sight and vr parallel to the line of sight (Fisher et al. 1994). It can be demonstrated that, if 
the ACF is expressed as ^(r) = (r/ro)"'*', then 

w{r,) = A{^)ryp-\ (1) 



^ImCrab in the 15-55 keV band corresponds to 1.27x10 ^'^erg cm ^ s 
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where ^(7) = r(l/2)r[(7 - l)/2]/r(7/2) (Peebles 1980). 

The ACF has been estimated by using the minimum variance estimator described by Landy 

& Szalay (1993): 

, DD-2DR + RR 

arp,rr) = — (2) 

where DD, DR and RR are the normahzed number of data-data, data-random, and random- 
random source pairs, respectively. Equation |2] indicates that an accurate estimate of the 
distribution function of the random samples is crucial in order to obtain a reliable estimate of 
^(rp, vr). Several observational biases must be taken into account when generating a random 
sample of objects in a flux limited survey. In particular, in order to reproduce the selection 
function of the survey, one has to carefully reproduce the space and flux distributions of the 
sources, since the sensitivity of the survey in not homogeneous on the sky. Simulated AGN 
were randomly placed on the survey area. In order to reproduce the flux distribution of the 
real sample, we followed the method described in MuUis et al. (2004). The cumulative AGN 
logN-logS source count distribution, in the 15-55 keV band, can be described by a power 
law, S = kS'"", with a ~ 1.55 (Ajello et al. 2008) Therefore, the differential probability 
scales as S'^*-"^^-'. Using a transformation method the random flux above a certain X-ray flux 
Slim is distributed as S = Sumi^ ~ 'P)~ 1 where p is a random number uniformly distributed 
between and 1 and S'/jm,=7.6xl0~^^ erg cm~^ s~^. All random AGN with a flux lower than 
the flux limit map at the source position were excluded. Redshift were randomly drawn 
from the smoothing of the real redshift distribution with a gaussian kernel with o"2=0.3. An 
important choice for obtaining a reliable estimate of w(rp), is to set Timax in the calculation 
of the integral above. One should avoid values of T^max too large since they would add noise 
to the estimate of w{rp). If, instead, iTmax is too small one could not recover all the signal. 
We have calculated w{r^) by varying r^max and found that the result converges at i^max ~60 
Mpc/h. 

Errors on w(rp) were computed with a bootstrap resampling technique with 100 realizations. 
It is worth noting that in the literature, several methods are adopted for errors estimates in 
two-point statistics, and no one has been proved to be the most precise. However, it is known 
that Poisson estimators generally underestimate the variance because they do consider that 
points in ACF are not statistically independent. Jackknife resmpling method, where one 
divides the survey area in many sub flelds and iteratively re-computes correlation functions 
by excluding one sub-fleld at a time, generally gives a good estimates of errors. But it 
requires that sufficient number of almost statistically independent sub-fields. This is not the 
case for our small sample. For these reasons we used the bootstrap resampling for the error 
estimate which, in our case, are comparable with the Poisson errors. 
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Results 



We show in Fig. [2] the projected ACF measured on the whole AGN sample of the 
survey. The ACF has been evaluated in the projected separation range ~0.2 Mpc/h < Vp < 
200 Mpc/h and has been plotted in form of w{rp)/rp in order to reproduce the slope of 
^(r) (see above). The bin size for computing w{rp) has been set to Alog{rp, 7c)=0.15. We 
obtained an estimate of w{rp)/rp with a significance of the order 4a-5a. In order to evaluate 
the power of the clustering signal we fitted w{rp) with x^ minimization technique by using 
EqU] as a model with tq and 7 as free parameters. The correction due to the integral 
constraint (Peebles 1980) is estimated to be much smaller than the statistical uncertainties 
in our sample, and thus has not been made. As a result we obtained ro=5.56lQ;43 Mpc/h and 
7=1.64^Qg°^. The confidence contours of the fit are presented in Fig. [2j We also measured 
the ACF for different data subsamples. We first divided the sample according to the column 
density: we defined as Type II AGN (or absorbed) sources with logi^Nn) >22 cm~^ and as 
Type I AGN (or unabsorbed) sources if \og{NH) <22 cm~^. As a result we constructed a 
sample of 96 Type I AGN and one of 103 Type II AGN. For both samples we computed 
the ACF with the technique described above. We also split the sample into high and low 
luminosity subsamples. All the sources with Li5_55 >43.2 erg/s (i.e. the median luminosity 
of the whole sample, HL sources) were classified as high luminous, while the sources with 
L15-55 <43.2 erg/s (LL sources) as low luminosity sources. The results of the measurement 
of the ACF as a function of the source type and luminosity class are presented in Fig. [2] 
together with the fit confidence contours. Note that for the HL sample the fit parameters are 
poorly constrained because of the lack of close pairs in the sample. We also repeated the fit 
by freezing 7 to 1.7, and obtained consistent results (Tabled]). A summary of the fit results 
of all the samples used here is given in Table [TJ Type I AGN show a larger correlation with 
respect to that of type II AGN, the significance of this difference is of the order ~ 2.7cr — 3.3cr. 
HL AGN show a 1.7-4.6o" higher correlation length with respect to LL AGN. We also checked 
the correlation between ro and < Lx > of all the subsamples and found a linear correlation 
coefficient R=0.95, which corresponds to a ~ 2o" significant correlation. We can use the 
weighted mean dispersion of the results on the measurement of ro in our subsamples to 
estimate the impact of sample variance on our results under the assumption that this is the 
main reason of the fluctuations. Our estimates suggest that overall our results are affected 
by this effect at ~10% level. It is worth to note that our results are more significant than 
those obtained by e.g. MuUis et al. (2004), with a similar number of sources. This is because 
our sources are distributed in a much smaller volume than that sampled by the NEP survey 
and, by being on average less luminous, have an intrinsic higher space density resulting in a 
larger number of close source pairs. 
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In the linear theory of structure formation, the bias factor defines the relation between 
the autocorrelation function of large scale structure tracers and the underlying overall mat- 
ter distribution. In the case of X-ray selected AGN, we can define the following relation: 
^x{i",z) = bx{r,zY^rn{r,z), where ^x, ^m and hx are the autocorrelation function of AGN, 
of DM and the AGN bias factor, respectively. In order to compute the hx, we estimated the 
amplitude of the fluctuations of the AGN distribution in a sphere of 8 Mpc/h (also know as 
(Tg), by using Eq. 12 and 13 of Miyaji et al. (2007) and all the combinations of ro and 7 are 
reported in Tab. [TJ In order to derive the bias factor of the AGN in our samples we used 
bAGN = o"8,agn(-2)/o"8-D(2;), where D(z) is the growth factor. This quantity allows us to com- 
pare the observed AGN clustering to the underlying mass distribution from linear growth 
theory (Hamilton 2001). As a result we obtain for the whole sample byiG7v(z~0.04)=1.21^QQ7. 
We have repeated this calculation for all the samples listed in Table [T] and we report the 
corresponding values of the bias factor in Table U\ for all the possible best fit results. 
It is widely accepted that the clustering amplitude of DMH depends primarily on their mass 
(see e.g. Sheth et al. 2001). In this way, we can estimate the typical mass of the DMH 
in which the population of AGN reside, under the assumption that the typical mass of the 
host halo is the only variable that causes AGN biasing. We have then computed the ex- 
pected large-scale bias factor for different dark matter halo masses by using the prescription 
of Tinker et al. (2005). The required ratio of the critical overdensity to the rms fluctuation 
on a given size and mass is calculated by z/ = Scr/cr^M, z), for our purposes we assumed 
6cr ~ 1-69 and compute cr8(M, z) and therefore b{M,z) using Eq. A8, A9, and AlO in van 
den Bosch (2002). The typical DMH mass that hosts an AGN has been estimated to be 
log{M]:)MH)=^'^-^5-o'i3 h~^M/M0. This is consistent with similar measurements in the local 
Universe of Krumpe et al. (2010), Grazian et al. (2004) and Akilas et al. (2000). We have 
computed the typical mass of the DMH for all the subsamples listed in Tableland reported 
for simplicity in Table [2l Following Martini & Weinberg(2001), by knowing the AGN and 
DMH halo space density at a given luminosity and mass {tiagn^iT'Dmh), one can estimate 



Table 1: Summary of the results. 
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"Number of sources in the Sample. 

''ro obtained by freezing 7=1.7 in the fit. 
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Table 2: Bias factor and Mass of the Dark matter halos hosts of the AGN in the samples. 



Sample 
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logiMBH)") 


logi^r 


M*^ 






log(M/ h^^Mg) 


Gyr 


log{M/MQ) 




10^°/Mq 
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, ,-,+U.OB 


,0 ,^+o.ua 


0.68 


8.51 


-1.96 


18.2 
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2.0ltVll 


10 Q/+C1.15 


4.99 


8.79 


-2.02 


31.6 


Type II 


1.08lS-fg 


l2.92+°;ii 


1.32 


7.96 


-1.85 


6.38 


HL 


9 90 + 0.95 


14.08±g;f„ 


3.91 


9.28 


-2.12 


80.5 


LL 


o.80_o;ig 


j^gg+a34 


0.24 


7.43 


-1.68 


2.28 



"AGN bias factor 

^Mass of the typical Dark matter halo hosting an AGN 

^AGN duty ciycle in Gyr 

'^Black hole mass 

'^Eddington ratio 

•^Stellar mass of the bulge 



the duty cycle of the AGN, TACNiz) = '^^^^^(jf^'^niz)- Where th{z) is the Hubble time 
at a given redshifto- For the whole sample at z~0, udmh ~6.7x10~^ Mpc"'^ (Hamana et 
al. 2002) and uagn ~3.4x10~^ Mpc~^ (Sazonov et al. 2007) which leads to an estimate of 
TAGNiz = 0) ~0.68 Gyr. To fully characterize our sample, we derived the average properties 
of the active BHs and their host galaxies. By using the bolometric correction prescribed 
by Hopkins et al. (2007) we estimated from < Lx >, < LboI > and < Lb > (B band 
luminosity). L^ is related to the black holes mass and the stellar mass of the host galaxy 
via scaling relations (Marconi & Hunt 2003). From < Mbh > we derived < LEdd > and the 
Eddington rate e=< LboI > / < LEdd > (see Table [2] for a summary). We point out that the 
estimates obtained above have several limitations which mostly arise from the uncertainties 
on scaling relations and from the broad range of luminosities sampled here. We therefore 
consider these values as estimate for the "average" local AGN population. 



5. Summary and discussion 

In this letter we report on the measurement of clustering of 199 AGN in the local 
Universe using the Swift/BAT all-sky survey sample. This result gives, for the first time, 
an unbiased picture of the z=0 DMH-galaxy-AGN coexistence/evolution. We obtained a 
correlation length ro=5.56lo;43 Mpc/h and 7=1.641q;q8. We measured the ACF for Type 



This is an upper limit obtained by assuming that the lifetime of the DM halo is of the order oi th{z) 



I and Type II AGN and found a significant difference in their correlation lengths. We 
have measured a marginally significant higher ro for high luminosity AGN than for the low 
luminosity ones. We propose that the observed difference in Type I vs. Type II clustering 
is driven by the intrinsic higher < Lx > of Type I AGN as we show a marginal evidence 
of a correlation between ro and Lx- We estimated the typical mass of the DMH hosting an 
AGN of the order log(MMD_ff)~13.15 h~^M/M0. In Fig. [3] we show the bias-redshift plane 
results from AGN and galaxy surveys (references in the figure). In the same plot we show 
the expected evolution of different DMH masses. We compared only bias values of studies 
that rely on the real space correlation function ^(r) (values of (Tsagn/gal from Krumpe et 
al. 2010). This approach allows us to compare all different clustering studies on a common 
basis. 

The majority of the X-ray surveys agree with a picture where AGN are typically hosted 
in DM halos with mass of the order of 12.5 h'^M/Mo < log(MA/D//)<13.5 h"^M/M0 which 
is the mass of moderately poor group. Optically selected AGN instead reside in lower density 
environment and of the order of the log(MM_DH)~12.5 h~^M/MQ. Another interesting aspect 
is that X-ray selected AGN samples (including ours) cluster similarly to red galaxies and 
that LL AGN or type II AGN are found typically in less massive environments. On the 
contrary HL AGN and Type I AGN are hosted in massive galaxies in massive DM halos 
(clusters). 

We estimated that Swift-BAT AGN are powered by black holes with a typical mass 
log(MBi//M0) ~8.5, accreting at very low Eddington ratio (i.e. ~0.01 LEdd) and that they 
are hosted by massive galaxies with mass of the order ~ 2 x 10^^ M/Mq. These properties, 
except e, scale with < Lx > and Type (HL and Type I are hosted in higher mass DMH in 
more massive galaxies with bigger black holes). The upper limits on the duty cycle suggest 
that these AGN are shining since at least for 0.2-1.2 Gyr|j. 

We then tested the AGN merger-driven triggering paradigm by comparing the theoret- 
ical predictions for AGN clustering of the model of Marulli et al. (2008) and Bonoli et al. 
(2009) with our results on the whole sample. The theoretical model is based on the assump- 
tion that AGN activity is triggered by galaxy mergers and the lightcurve associated to each 
accretion event is described by an Eddington-limited phase followed by a quiescent phase 
modeled after Hopkins et al. (2005). Applied to the Millennium Simulation (Springel et al. 
2005), such model has been shown to be successful in reproducing the main properties of the 
black hole and AGN populations (Marulli et al. 2008) and the clustering of optical quasars 
(Bonoli et al. 2009). Using this model, we computed the expected correlation length of a 



'HL and type I estimate may be wrong because of the relatively young age of 10 (M/Mq) DMH 
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sample of simulated AGN at z~0 with intrinsic Lboi luminosities similar to the ones of our 
observed AGN. The model predicts a clustering length ro=5.68±0.08 which is in agreement 
within Icr with our measurement. 

By merging the observational evidences and the model predictions, a plausible scenario 
for the history of local AGN is the following: 

• Swift- BAT AGN switched on about 0.7 Gyr Ago after a galaxy merger event. 

• They shine in an Eddington-limited regime for the first part of their lives where they 
gain most of their mass. 

• In the second phase of their lives (i.e. after 0.2-0.5 Gyr) they start to accrete with lower 
and lower efficiency. Their luminosity drops because of the decreased gas reservoirs. 

• At z~0 they have grown to ~10^~^ M/Mq SMBHs, shining as moderately low-luminosity 
AGN at low accretion rates. 
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Fig. 2. — Upper Left Panel : w{rp)/rp for the whole sample, the solid line represents the 
best fit power-law. Upper Right Panel : The ACF of type I AGN (blue) and Type II (red). 
Bottom Left Panel : The ACF of HL [blue) and LL {red) AGN. Bottom Right : Confidence 
contours of the two parameter fit to ^(r) for the whole sample (black), Type I AGN (green), 
type II AGN (pink), HL AGN (blue) and LL (red). 
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Fig. 3. — Bias parameter 6agn as a function of redshift for studies that are based on the real 
space correlation function. Yellow (green) symbols represent high (low) luminosity subsam- 
ples of the corresponding study. The blue points for galaxy clustering measurements belong 
to blue galaxies, while red points represent the red galaxy sample. The dotted lines give the 
expected b{z) of typical dark matter halo masses based on Sheth et al. (2001) and van den 
Bosch (2002). The masses are given in log Mdmh in units of h~^ Mq. 
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